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Abstract 

This paper presents a study of vertex tracking with a beam hodoscope consisting 
of three layers of monohthic pixel sensors in SOI technology on high-resistivity 
substrate. We study the track extrapolation accuracy, two-track separation and 
vertex reconstruction accuracy in vr^ Cu interactions with 150 and 300 GeV/c pions 
at the CERN SPS. Results are discussed in the context of vertex tracking at future 
lepton colliders. 
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1 Introduction 



Vertex tracking is one of the capabilities most crucial to a detector at future 
lepton colliders and flavour factories. Much of the anticipated physics program 
at a high energy lepton collider, such as the ILC, CLIC or a muon collider 
relies on the ability of efficiently discriminating heavy {t, b, c and r) from 
light fermions [112] • If the Higgs boson exists at a mass of ^ 125 GeV, as 
possibly indicated by the recent preliminary LHC [31H] and Tevatron [5] re- 
sults, the precise determination of its couplings to fermions will be essential 
to test whether the Brout-Englert-Higgs mechanism of electro-weak symmetry 
breaking is also responsible for fermion mass generation. Experiments at lower 
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energy facilities, such as the Super KEKB [6] and SuperB |7] 6-factories, re- 
quire high resolution vertex tracking for time-dependent measurements, which 
are crucial to the physics program in heavy flavour physics [8]. Because the 
momentum of particles from B meson decays at the T(4S) is low and the sen- 
sitivity to short-lived hadrons produced at high energy colliders must extend 
to most, if not all, of their charged decay products, excellent track extrapola- 
tion accuracy at low momentum is of paramount importance. This motivates 
the development of thin Si sensors. Monolithic pixel sensors, such as CMOS 
active pixel sensors and DEPFET sensors, have emerged from the R&D effort 
for future e'^e~ colliders as the best suited solution for high-resolution ver- 
tex tracking with a minimal material budget in a low to moderate radiation 
environment PllOllllj . More recently, monolithic pixel sensors on high resis- 
tivity substrate have demonstrated superior performance owing to the larger 
amount of charge collected, ensuring high detection efficiency, and the faster 
collection time, reducing also the effect of radiation damage [T2II3]- 

Among the technologies offering a CMOS process over a high resistivity sen- 
sitive volume, Silicon-On-Insulator (SOI), pioneered a decade ago [H], has 
been extensively tested in the last six years, within a collaborative R&D ef- 
fort developed in partnership with KEK and Lapis Semiconductor Co. Ltd. 
(formerly OKI Semiconductor). With the mitigation of the back-gating effect 
using a buried p-well (BPW) to protect the in-pixel CMOS circuitry |I5], SOI 
pixel sensor prototypes have demonstrated high detection efficiency and sin- 
gle point resolution at the micron level [16] , also after thinning of the handle 
wafer [T7] . 

In this paper we illustrate the vertex tracking performance of a hodoscope 
composed of SOI pixel sensors based on the study of inelastic interaction of 
high energy pions in a beam test (SOIPIX) at the CERN SPS in Fall 2011. 
We present the results obtained for tracking resolution, two-track separation 
and vertex resolution, compare our data to simulation and discuss the relation 
to requirements for future colliders. 



2 Experimental Setup and Data Analysis 

The SOIPIX beam test experiment has been carried out on the SPS beam-line 
H4 in the CERN North Area with a beam hodoscope made of three layers of 
"SOImager-2" prototype chips, designed at LBNL and produced in 0.2 
Lapis technology (formerly OKI) on n-type SOI wafers with a resistivity of 
the handle wafer of ~ 700 Q-cm. The sensor sensitive area is a 3.5x3.5 mm^ 
matrix of 256x256 pixels arrayed on a 13.75 /xm pitch, read out through four 
parallel arrays of 64 columns each [16] . The sensor detection performances for 
high energy particles and soft X-rays have already been presented in details 
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in jl6fl7fl8j . In this setup, a single sensor, thinned to 75 /im (singlet), is 
located upstream from a Cu target to define the position of impact of the 
incoming pion. A pair of 260 /xm-thick sensors spaced by 9.4 mm (doublet) 
reconstructs the trajectories of the charged particles emerging from the target. 
This spacing is comparable to that of the two layers of the PXD pixel detector 
in BELLE-II and to the distance between the first and second layers of a LC 
vertex tracker. The 3 mm-thick Cu target is inserted between the detector 
singlet and the doublet, 15 mm upstream from the first layer of the SOI pixel 
doublet. The extrapolation distance from the first doublet layer and the target 
corresponds to the design distance between the interaction point and LayerO at 
SuperB and the first layer of the vertex tracker of an e'^e~ linear collider at ^/s 
= 250 - 500 GeV, while for a multi-TeV e~^e~ collider, such as CLIC, the beam 
induced background requires a stay-clear radius of ~25-30 mm. The use of a 
thin sensor for the singlet minimises the interaction effects upstream from the 
target. Data have been collected with 7r~s of 150 and 300 GeV/c momentum. 
Runs with 300 GeV/c pions have also been taken without the target, for 
alignment and calibration purposes. We operate the sensors at a depletion 
voltage, Vd, of 50 V, corresponding to a depleted thickness of ~100 /im in the 
doublet sensors and ~60 /xm in the thin singlet sensor. 

The data acquisition system consists of a custom analog board pigtailed to a 
commercial FPGA development board, equipped with a Xilinx Virtex-5 FPGA 
used as control board [I9]. The sensor analog outputs are fed to 100 MS/s 
14-bit ADCs through independent analog differential inputs. Digitised data 
are formatted and transferred to the DAQ computer via a USB-2.0 link at a 
rate of 25 Mbytes/s. Measurements are performed with the chip clocked at 
12.5 MHz, corresponding to an 80 ns read-out time per pixel. The readout is 
synchronised with the SPS extraction pulse, such that a sequence of 960 frames 
from each of the three sensors are read-out during the 9.6 s-long SPS spill. 
Detectors are cooled using forced airflow and the temperature near the sensor 
surface is in the range (28±1)°C during data taking. Data sparsiflcation and 
zero suppression is performed on-line using a custom Root-based [20] program. 
Sensors are scanned for seed pixels with signal exceeding a preset threshold 
in noise units. For each seed, the 7x7 pixel matrix centred around the seed 
position is selected and stored on flle. The pixel pedestal and noise values are 
updated at the end of each SPS spill, in order to follow possible drifts of their 
baselines throughout the data taking. Data are stored in Root format and 
subsequently converted into Icio format \2T\ for offline analysis. 

The data analysis is based on a set of custom processors in the Marl in recon- 
struction framework [22] to perform cluster centre-of-gravity reconstruction, 
pattern recognition, track and vertex fltting [23]. Clusters are reconstructed 
applying a double threshold method on the matrix of pixels around a selected 
candidate cluster seed. Clusters are requested to have a seed pixel with a 
signal-to-noise ratio, S/N, of at least 7.0 and the neighbouring pixels with a 
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S/N in excess of 5.0. Clusters consisting of a single pixel are discarded. The 
cluster position is computed using the centre of gravity of the measured pulse 
height. The observed average and most probable signal-to-noise ratio for clus- 
ters associated to reconstructed particle tracks in the pixel doublet is 47.4 and 

43.0. The most probable value of the signal-to-noise ratio for the seed pixel is 

45.1. This can be compared to 45.2±0.4 predicted by sensor simulation [T7] . 
The detector planes have been optically surveyed after installation on the 
beam-line. Their measured positions are used as starting point for the offline 
alignment procedure with particle tracks, performed using the millipede-2 
program |24] . 

A total of 201000 and 220000 events, having at least one hit per layer in the 
doublet and only one hit in the singlet, have been recorded in the September 
2011 data taking with the target installed and 150 and 300 GeV/c beam, re- 
spectively. An additional ~20000 calibration events have been collected with- 
out the target for alignment. 

Simulation of interactions and energy release by charge particles in the detec- 
tors is performed using the Geaiit-4 simulation toolkit |25] with the FTFP_BERT 
physics list [26], which implements high energy inelastic scattering of hadrons 
by nuclei using the FRITIOF model [27]. The yield of short- and long-lived 
particles has been studied using Pythia 6.125 [2B]- Charge collection and 
signal generation in the pixels is simulated using a custom Marl in processor, 
PixelSim [22], where the pixel S/N is tuned to the values obtained in data for 
each chip. Simulated events are then analysed through the same reconstruction 
chain as real data. 



2. 1 Track Reconstruction 



The efficiency of the SOI sensors is determined by using particle tracks recon- 
structed on the other two layers and extrapolated to the layer under study, in 
the events taken without the target. 

The Geant-4 + PixelSim simulation predicts a sensor efficiency of O.OOSlg;!}?! 
for the doublet. We estimated the efficiency of the sensors in the doublet in 
an earlier beam test performed in 2010, when the same SOI sensors were used 
both in for doublet and the singlet. We use the number of tracks reconstructed 
from hits on the second layer of the doublet and on the singlet which had an as- 
sociated hit on the first layer of the doublet and found a value of O.OSIq oI [TB] . 
The single point resolution of the sensors is measured to be (1.12±0.03) /xm 
for the doublet [16] and (1.7 ±0.50) /im for the singlet [17], where simulation 
predicts (1.07±0.04) /im and (1.63±0.05) /xm, respectively. The difference is 
due to the lower amount of charge collected, and thus of signal-to-noise ratio. 
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in the thin sensor. 



Inelastic interaction events are reconstructed from tracks formed with the 
space points in the two doublet layers. We select events with only one recon- 
structed hit in the singlet layer and at least one accepted hit in each of the 
doublet layers to perform the track reconstruction. This pre-selection removes 
empty events and most of those with more than one primary pion recorded in 
the same frame. 

We reconstruct the trajectory of the incoming pion by projecting the posi- 
tion of the hit reconstructed in the thin singlet. Since the r.m.s. divergence of 
the beam is 1.5x10"^ rad in the vertical and 7.5x10"^ rad in the horizontal 
coordinate, as measured in events taken without the target, the uncertainty 
on the vertical extrapolation of the pion position of impact on the target is 
^3 /im r.m.s. Secondary candidate tracks are reconstructed in the doublet 
by pairing the hits. These are extrapolated to the singlet and their impact 
parameter w.r.t. the hit recorded on the singlet sensor computed. Its dis- 
tribution is peaked at the value of the extrapolation resolution for primary 
hadrons, which have not interacted in the target, while it has larger values 
for secondary particles produced in interactions in the target, as well as for 
primary pions experiencing large scattering, as shown in Figure [H We resolve 
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Fig. 1. Impact parameter resolution of tracks reconstructed in the SOI doublet and 
extrapolated to the singlet hit (in mm). The distribution of tracks from inelastic 
events are shown by the filled histogram. 

the combinatorial in the pairing of the doublet hits by choosing the hit pair 
minimising the impact parameter, provided that the longitudinal position of 
the intercept of the corresponding track with the estimated primary hadron 
track is contained in the fiducial region between the singlet layer and the first 
doublet layer. In simulation, this strategy correctly assigns the hit pairing for 
0.96 of the particle tracks emerging from an interaction and contained within 
the geometrical acceptance of the detector. After this first vertex is 

reconstructed for events with two or more tracks. Then, candidate tracks are 
formed by pairing the unassociated hits in the doublet and tested for compat- 
ibility with the interaction vertex by computing their impact parameter w.r.t. 



5 



its position. Candidate tracks with impact parameter smaller than 40 jjm are 
kept. Again, the combinatorial is solved by using the hit pair minimising this 
impact parameter. 



2.2 Vertex Reconstruction 

The interaction vertex for inelastic interaction events is reconstructed with the 
candidate tracks fitted in the doublet. Several algorithms have been designed 
to reconstruct the interaction vertex from the trajectories of charged particles. 
They are typically implemented in the context of collider experiments with a 
homogeneous magnetic field using a helix track model. Due to the geometry 
of the beam test setup and the lack of magnetic field, these implementations 
have numerical problems in reconstructing the interaction vertex for our ex- 
perimental setup. In absence of a magnetic field, the track model is a simple 
straight line with four parameters. While the track model is simple, a specific 
implementation of the vertex fit procedure is required. Due to the low multi- 
plicity in the interaction events, the vertex position {xy,yy,z^) and the track 
parameters Pi = {x'^^y'^ = {{dx/dz)i,{dy/dz)i) for particle i can be simply 
derived from the minimisation of the sum of contributions of the trajectory 
to a global fit. The measurement of the trajectory for particle i on detector 
i at position is given by its Cartesian coordinates {xj,yl) and the corre- 
sponding reconstruction errors and ay. The xl contribution of particle i is 
then defined as 



Monte Carlo simulation studies show that this approach avoids the numer- 
ical shortcomings of other implementations of vertex fits. The choice of an 
adequate track model is essential to reach high efficiency and accuracy in the 
fitting procedure. In this implementation a global fit has been used for sim- 
plicity, but the straight line track model could also be implemented in more 
complex vertex reconstruction programs. The vertex fit is performed itera- 
tively. If the vertex probability is smaller than 10~^, the track giving the 
largest contribution is removed and the vertex recomputed. The hits used 
for reconstructing these tracks are not further used. A reconstructed event is 
shown in Figure |5J Using fully simulated and reconstructed events, we estimate 
the vertex reconstruction efficiency for inelastic interaction events with at least 
two charged particles within the detector acceptance to be 0.94±0.05. Further, 
the reconstructed multiplicity of tracks associated to the vertex agrees with 
the generated charged multiplicity in the detector acceptance within ±0.04. 



j 




(1) 
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Fig. 2. A 300 GeV/c vr" inelastic interaction event with five charged reconstructed 
particles. The beam comes from the left and the pion position is first determined in 
the singlet layer, shown on the left. After the interaction in the target, in the centre 
of the picture, charged particle tracks are reconstructed on the two doublet layers 
shown on the right. 

3 Results 



The number of candidate inelastic interaction events, reconstructed in data by 
the procedure discussed above, is 1105 and 2118 for 150 and 300 GeV/c pions, 
respectively, corresponding to 0.0055±0.0002 and 0.0096±0.0006 of the events 
preselected. The observed fractions agree quite well with simulation, which 
predicts 0.0065±0.0003 and 0.0101±0.0003, confirming the efficiency indicated 
by simulation. For the subsequent analysis we apply additional quality cuts on 
the reconstructed events requiring the vertex fit probability to exceed 10~^ 
and the uncertainty on the vertex longitudinal position to be smaller than 
1.5 mm. This yields a total of 1048 selected interaction events at 150 GeV/c 
and 1862 at 300 GeV/c. 



3.1 Tracking Resolution 



The track extrapolation accuracy at the position of the Cu target, located 
15 mm upstream from the first layer of the SOI pixel doublet is 3.7 /im for 
300 GeV/c and 8.9 /im for 1 GeV/c particles according to simulation, as 
shown in Figure [31 These values are quite comparable to those expected for 
the vertex tracker proposed for a linear collider, i.e. ~2 /xm p^hTe^ — ' 
where 9 is the track polar angle, for ^/s = 250 - 500 GeV with innermost 
pixel layer at ~15 mm radius, and ~4 /im ^^ ^^iit^e^ — 5 ^r y/s = 3 TeV with 
innermost layer at ~30 mm radius. The impact parameter of tracks w.r.t. the 
reconstructed vertex is computed for inelastic interaction events with three or 
more tracks associated to the vertex. Tracks are iteratively removed from the 
vertex and its position recomputed to avoid the bias due to the contribution 
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Fig. 3. Track extrapolation resolution on the target using the SOI doublet. Left: 
the simulation result is shown by the dots with error bars while the fitted curve is 
3.7 /im ^.0 fimGeY — ^ Right: impact parameter of tracks w.r.t. the interaction ver- 
tex. Tracks are removed from the vertex before the impact parameter is computed. 
Points with error bars represent the data and the line the simulation. 

of the track under study. The distance between the track and the interaction 
vertex is computed at the position of closest approach. The most probable 
value of the impact parameter resolution is 12.7 /xm, which is the convolution 
of the track extrapolation resolution, the vertex position resolution (~9.0 fim, 
see section 3.4) and the multiple scattering in the target. Data and simulation 
agree well, as shown in Figure [31 

3.2 Two-track Separation 

Two-track separation is an important figure of merit for vertex tracking and 
the optimisation of the pixel sensor geometry. In fact, a given target value 
of the single point resolution can be achieved with different pixel sizes, P, 
depending on the S/N ratio and the charge spread on neighbouring pixels. In 
the case of the sensor used here, with P = 13.75 /im, a seed most probable S/N 
value of 45 and ~0.35 of the total cluster charge collected on the surrounding 
pixels, we obtain a 1.1 /xm resolution for isolated particle tracks. The same 
resolution could be achieved also with a pitch P = 25 /xm, for a S/N ~ 70. 
On the other hand, two-track separation, which scales oc P, and occupancy, 
which scales oc P^, put emphasis on small pixel size and limited charge spread 
to minimise cluster merging. In inelastic interaction events we observe a hit 
density on the first doublet layer of 2.5 hits mm~^. The two-track separation 
is studied on these events by computing the distance between the two closest 
reconstructed hits on the two layers of the doublet. Due to the large incoming 
pion momentum, the interaction products are highly boosted in the forward 
direction, closely resembling a highly coUimated jet at collider experiments. 
When they are detected on the first layer of the doublet, the minimum hit 
distance has a most probable value of 70 fim. In the core of simulated hadronic 
jets it is ~400 /im in e~^e~ — )■ H^Z^ — )■ bbqq events at y/s = 500 GeV and 
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Fig. 4. Two track separation characterisation. Left: minimum distance between a 
hit associated to a reconstructed track and the closest hit on the first layer of the 
SOI doublet. Right: number of pixels along rows in clusters associated to tracks. 

~240 /iin in e'^e^ — )• H^A'^ — )• bbbb events at ^/s = 3 TeV. In our data, we 
observe a distinct cut-off of the distribution at a minimum hit distance value 
of ~50 fim, corresponding to the smallest two-track separation resolved in 
the detector. This is comparable to the average size of clusters associated to 
reconstructed tracks along rows and columns, which is 2.67±0.6 pixels, and to 
its most probable value, 3 pixels corresponding to 42 /im (see Figure H]). Below 
this distance, the charge deposited on the pixels by distinct tracks is merged 
into a single cluster hit. We study the hit merging by selecting events with a 




Cluster Pulse Height (ADC Counts) 

Fig. 5. Two track separation and hit merging. Left: display of event with two tracks 
(labelled 2 and 3) giving a merged single hit with cluster pulse height of 1200 ADC 
counts on the first layer of the SOI doublet. Right: Pulse height for clusters on the 
first layer of the doublet in events with a pseudo-track made with an unassociated 
hit in the second layer (unfilled histogram). The grey filled histogram shows the 
pulse height distribution for clusters which are within 50 //m from the pseudo-track 
extrapolation and are already associated to another track. The continuous line shows 
the fit of a Landau function. 

reconstructed interaction vertex and at least one hit in the second SOI doublet 
layer which is not associated to a track. In these events, a "pseudo-track" is 
built by fitting a straight line from the unassociated hit to the reconstructed 
vertex. This pseudo-track is extrapolated to the first layer of the doublet and 
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the hit closest to the extrapolation point is selected. Figure |5] shows one of 
such events and the cluster pulse height of all these hits with highlighted those 
having a distance from the extrapolated point smaller than 50 /im, i.e. three 
and half times the pixel pitch, and which are already associated to another 
reconstructed track. This category is characterised by a most probable value 
of the cluster charge of (772±80) ADC counts, to be compared to 395±15 for 
the distribution of all the selected hits, i.e. larger by a factor of 1.95±0.02. 
This shows that tracks impacting on the detector at a distance closer than 
3-4 pixels produce merged hit clusters, characterised by a most probable value 
of the collected charge which is approximately double compared to that of 
isolated clusters. This characteristic can be used, in principle, for identifying 
them. 

It is useful to compare these results to the anticipated track density and two- 
track distance at high energy lepton colliders. These are driven by tracks in 
collimated hadronic jets and the amount of machine-induced backgrounds, 
mostly incoherent pairs and 77 — )■ hadrons events, integrated in a read-out 
cycle. In the process e"'"e~ — )■ H'^A'^ — )■ bbbb at a/s = 3 TeV, the fraction of 
particles from a b hadron decays reaching the first layer of the Vertex Tracker, 
at a radius of 30 mm, and within less than 50 /^m from another charged particle 
track is 0.03. The probability to have a charged particle from beam-induced 
background closer than 50 /zm is 0.02 assuming a detector integration time of 
25 ns, corresponding to 50 bunch crossings. 



3.3 Interaction Multiplicity 

The charged multiplicity of inelastic interactions is measured from the number 
of particles associated to a reconstructed vertex, as discussed above. In order 
to compare with existing data, we include one-prong interaction events. These 
are events having a single reconstructed track in the SOI doublet, which has a 
longitudinal intercept with the primary pion trajectory in the range 16 < z < 
23 mm, i.e. compatible with the target position, and an impact parameter 
w.r.t. the position of the hit in the SOI singlet larger than 100 fim. We obtain 
a muhiphcity of 2.16±0.18 and 2.57±0.22 for 150 and 300 GeV/c vr" beam, 
respectively (see Figure |6]). The measured multiplicity scales logarithmically 
with the beam energy, as expected [30]. We compare the multiplicity measured 
in our data as a function of s to a fit to published data scaled to the rapidity 
range covered by our detector. Following [3lj, we use a function of the form 
< n > = a + b ln(s/so), where the parameters a and b are fitted on results 
reported for vr Cu collisions at beam energies from 50 to 200 GeV/c [32]. We 
get a=2.32, 6=0.84 for so=392 GeV^. Our results agree with this curve, as 
shown in Figured 
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Fig. 6. Reconstructed charged multiplicity of inelastic interaction events for 150 
(left) and 300 GeV/c 7r~s (right). Data are shown by points with error bars and 
simulated and reconstructed events by the continuous line. 
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Fig. 7. Average charged multiplicity of inelastic interaction events as a function of 
s. Data are shown by points with error bars and the continuous line represents the 
fit to earlier data rescaled to the same rapidity interval. The error bar associated 
to the line represents the uncertainty of the points used for the fit. 

3.4 Vertex Resolution 



The resolution on the vertex position reconstruction in the plane transverse to 
the beam is obtained by comparing the reconstructed vertex position to that 
of the incoming pion detected in the thin sensor upstream from the target. 
We measure a Gaussian width of (8.9±0.2) fim. By subtracting in quadra- 
ture the incident pion extrapolation resolution of 3 /im (see Section 2.1), we 
obtain (8.4±0.2) /im, which accounts for the vertex resolution, the multiple 
scattering of the incident pion in the target before the interaction and that 
of the emerging interaction products. The estimated average uncertainty from 
the vertex fit error is 6.2 /im and the difference between the generated and 
reconstructed position on simulated events has a Gaussian width of 6.5 /im. 

The resolution on the longitudinal vertex position is extracted first from the 
distribution of the longitudinal position of the vertex reconstructed in data 
and simulation, shown in Figure M Since the target has sharp edges, in ab- 
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scncc of resolution effects this distribution is a box function. Resolution ef- 
fects introduce a smearing of the box edges, which can be observed in both 
data and the simulated and reconstructed interaction events. These vary de- 
pending on the momentum and multiplicity of the detected particles emerg- 
ing from the interaction and on the longitudinal position of the interaction 
in the target. We perform a multi-parameter fit to the distribution of 
the reconstructed vertex longitudinal position. As a consistency test, we ob- 
tain a measurement of the target thickness of (3.15±0.10) mm in data and 
(3.07±0.05) mm in simulation, where a 3.00 mm thickness was generated. The 
average longitudinal vertex position resolution from the fit at 300 GeV/c is 
(503±60) fim and (532±39) /xm, again in data and simulation respectively. It 
becomes (615±83) fim and (590±35) fim for data and simulation, respectively, 
at 150 GeV/c where the multiplicity is lower and the multiple scattering is 
more important. 

We observe that the distribution of the vertex position away from the target 
is asymmetric, as we reconstruct more vertices downstream from the target 
than upstream. Since resolution and pattern recognition failures should give 
spurious vertices on either sides of the target position with equal probability, 
we interpret this as due to secondary vertices of long-lived particles produced 
in the primary interaction and decaying away from the target. Simulation 
shows a similar effect, mostly due to Kg tt+tt" decays, at a comparable 
rate as in our data. 

In data, the vertex fit error gives a most probable value of (335±7) /im and an 
average value of (490±6) fira, for the longitudinal vertex position resolution at 
300 GeV/c. The difference in the position of the generated and reconstructed 
vertices in simulated events has a Gaussian width of (290±10) /xm and non- 
Gaussian tails extending further, with the contribution of multiple scattering 
in the target to the Gaussian resolution estimated to be 60 /xm, at 300 GeV/c. 
Given the small angular coverage of our hodoscope and the large boost of 




Fig. 8. Distribution of the longitudinal position of fitted vertices in inelastic interac- 
tions of 300 GeV/c 7r~s. The points with error bars show the distribution for data 
and the continuous line for simulated and reconstructed events. 
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the interaction products, the resolution on the longitudinal vertex position 
is driven by the small opening angle of the detected interaction products. 
In 500 GeV e~^e~ collisions, simulated vertices from B mesons with energies 
in the range 100-150 GeV, re-weighted by the associated track multiplicity 
to match that of the reconstructed vertices in our data, give a resolution of 
(170±20) yum along the heavy meson line of flight [23]. On the other hand, 
at a multi-TeV collider the typical resolution for high energy jets is larger 
by a factor of two. Decay vertices in e~^e~ — )■ if°y4° — )■ bbbb events at ^/s = 
3 TeV, where the B mesons have energies in the range 200-700 GeV, have 
been reconstructed using the ZVTOP topological vertex algorithm [53]. The 
reconstructed secondary decay vertex has an average charged multiplicity of 
2.36 and the most probable and average values of the vertex position resolution 
along the jet axis are (330±5) /im and (530±3) /im, respectively. These values 
are very close to those obtained in the present study. 



4 Conclusions 

Tracking and vertexing performance of a small beam hodoscope made of SOI 
pixel sensors have been determined in a beam test with 150 and 300 GeV/c 
7r~s at the CERN SPS. The charged multiplicity of inelastic vr-Cu interac- 
tions has been determined to be 2.16±0.18 at 150 GeV/c and 2.57±0.22 at 
300 GeV/c, which agrees with earlier data and logarithmic scaling with the 
beam energy. The two-track separation of the sensor is ~50 /im, i.e. four pix- 
els, which agrees with the observed pixel multiplicity along columns and rows 
in reconstructed clusters. Merged hits can be identified by the large collected 
charge in the cluster. The vertex resolution is estimated to be 8.4 fim trans- 
verse and ^350-500 fim longitudinal. The results obtained are representative 
of the experimental conditions and reconstruction performance expected in 
highly boosted hadronic jets at TeV-class lepton colliders and demonstrate 
that pixels sensors with ~15 /im pitch on high-resistivity Si are extremely 
well suited for these applications. 
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